Influence of the photo-nuclear effect on electron-neutrino-induced electromagnetic 
cascades under the Landau-Pomeranchuk-Migdal regime in standard rock 
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The observation of earth skimming neutrinos has been proposed as a rather sensitive method to 
detect ultra-high energy (UHE) cosmic neutrinos. Energetic cosmic neutrinos can interact inside 
the rock and produce leptons via a charged current interaction. In the case of an incoming v e 
undergoing a charged current interaction, the produced UHE electron will induce an underground 
electromagnetic shower. At high energy (above 7.7 TeV in standard rock), such showers are subject 
to LPM (Landau, Pomeranchuk and Migdal) suppression of the radiative processes cross sections 
(bremsstrahlung and pair production). The consequence of this suppression is that showers are 
elongated. This effect will increase the detection probability of such events allowing deeper showers 
to emerge with detectable energies. On the other hand, the photonuclear processes which are 
usually neglected in electromagnetic showers with respect to radiative processes, turn out to become 
dominant in the LPM regime and will reduce the shower length. In this work, we have performed 
a complete Monte Carlo study of an underground shower induced by UHE electrons by taking into 
account both the LPM suppression and the photonuclear interaction. We will discuss the effects 
of both of these processes on the shower length and on the detectability of such events by ground 
arrays or fluorescence telescopes. We show that limits on neutrino fluxes that were obtained using 
simulations that were obviously neglecting photonuclear processes are overoptimistic and should be 
corrected. 



I. INTRODUCTION 

The search for ultra-high energy (UHE) (E> 10 17 eV) 
neutrinos requires a huge volume of target materials be- 
cause of their low interaction probability and low fluxes. 
The earth provides different target materials: some ex- 
periments are searching for UHE v interacting in sea wa- 
ter (for example ANTARES [1]), others in ice (IceCube 
[2-4] and ANITA [5, 6]) or in the atmosphere (HiRcs 
[7] and the Pierre Auger Observatory [8, 9]). One can 
also search for UHE v T interacting in the earth's crust 
in trying to detect r decay- induced showers [10]. Indeed, 
up-going v T can interact near enough to the earth's sur- 
face to produce a r whose range is long enough to es- 
cape the earth, decay and produce a shower in the atmo- 
sphere [11]. This shower can trigger a surface detector 
array such as the Pierre Auger Observatory Surface De- 
tector (SD) [12]. 

Another way of detecting up-going neutrinos is via the 
interaction of v e into the earth's crust. UHE v e will inter- 
act via charged current (CC) interaction with a nucleus 
of the medium and produce an electromagnetic and a 
hadronic shower. Since most of the incident neutrino en- 
ergy is transferred to the electron (^ 80% of E v ), the elec- 
tromagnetic shower will be very energetic. The hadronic 
shower produced by the neutrino CC interaction will not 
contribute significantly to the detection probability be- 
cause the energy transferred to it is negligible. In this 
paper we will discuss the development in rock of the v e 
induced electromagnetic showers. 
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An important effect plays a role in electromagnetic 
showers at high energy: the Landau, Pomeranchyk, and 
Migdal (LPM) effect, first predicted by Landau and 
Pomeranchuk [13] and derived in a modern quantum for- 
malism by Migdal [14, 15]. This effect induces an in- 
crease of the radiation length and hence the shower will 
be more penetrating and will have a higher probability 
to emerge in the atmosphere with large enough energy to 
be detected by a surface detector. 

Another effect needs to be taken into account which 
limits the LPM interaction length enhancement at high 
energy: the photonuclear interaction. The photonuclear 
effect consists of the interaction of a photon with a nu- 
cleus of the medium (via a direct or resolved process) 
which will produce a hadronic sub-shower. This effect, 
whose cross section increases logarithmically with energy 
of the photon, will increase the photon total cross section 
and therefore reduce its mean free path. On the other 
hand, the produced hadronic subshower will not prop- 
agate too far because of the small hadronic interaction 
length in rock. 

In this paper, we first start by presenting the LPM and 
the photonuclear cross sections and their energy depen- 
dence. We then discuss the consequences of these two 
effects on underground electromagnetic showers at UHE. 
We then illustrate this showing the results of a Monte 
Carlo study of the development of realistic showers in 
rock. 

We will conclude this work by studying the shower 
length and the probability that the emerging shower tail 
exceeds a given detection threshold. In this way we get 
a sense on their detectability by ground arrays or fluo- 
rescence detectors. We show that by neglecting photonu- 
clear processes one overestimates the detection sensitiv- 
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ity by a large factor and that earlier neutrino limits that 
were inferred in this manner are probably overoptimistic. 



II. ELECTRON AND PHOTON 
INTERACTIONS AT UHE 

Bcthe and Hcitler first calculated the cross sections of 
pair production and Bremsstrahlung. Some corrections 
were then added to these cross sections to take into ac- 
count screening effects. 

At UHE, the bremsstrahlung process takes place 
within a relatively large distance which is known as the 
photon formation length If and is given by If = h/q\\ 
where q\ | is the longitudinal momentum transferred to the 
target nucleus by an electron with energy E and mass m 
emitting a photon with energy k: qu = km 2 c 3 /2E(E—k), 
where c is the speed of light. If rises with electron energy. 
When If becomes larger than the interatomic spacing, 
then the amplitudes of multiple interactions with mul- 
tiple atoms add coherently and destructive interferences 
occur and suppress the bremsstrahlung amplitude (note 
that this was not taken into account in the original Bethe- 
Heitler theory). A similar reasoning applies to the pair 
production process [16]. The bremsstrahlung suppression 
is significant for the emitted photon with energy k from 
an electron with energy E when 

k , E{E k) 

»LPM 

with -Elpm given by 

m 4 A 



with E s = mc y/An/a = 21.2 McV, where a is the fine- 
structure constant. In standard rock (constant density 
p = 2.65 g.cm~ 3 , X = 27.6 g.cm -1 = 10.4 cm), one gets 
Elpm = 77 TeV so that in our case where we consider 
E 3> TeV, this effect cannot be neglected. The modifica- 
tion of the electron interaction length under LPM regime 
can be approximated by 

/ 2E 

X e = X@ \ I , 

V Elpm 

and for the photon by: 

Xl =X °\I' 50£ LPM ' 

For example, an electron with an energy of 10 18 eV will 
have its interaction length increased from 0.1 m to 22 m 
in standard rock because of the LPM suppression: the 
LPM effect will increase the shower length as mentioned 
in [17]. A purely electromagnetic shower under the LPM 
regime will develop over a much larger distance. 

Another process, playing a dominant role at high en- 
ergy, is the photonuclear effect. Photons can interact 



with nuclei of the medium and produce a hadronic sub- 
shower. At these energies, a photon interacts either di- 
rectly with a parton of the nucleus or may fluctuate into 
a quark antiquark state which will interact with the nu- 
cleus. Models based on different descriptions give fairly 
different cross sections for jp interactions. One can cite, 
for example, models such as in [18] assuming generalized 
vector meson dominance together with direct photon- 
quark interaction which takes into account resolved and 
direct photon interactions, model based on VDM (vec- 
tor meson dominance) with a Gribov- Glauber approach 
[19, 20] or the Donnachie-Landshoff model (exchange of 
pomeron pairs) in the Regge theory framework [21]. A 
comparison between Gribov- Glauber and GVDM may be 
found in [19]. 

In lower energy showers, this process is hidden behind 
the dominant radiative processes. In the LPM regime, 
this interaction becomes important as the cross section 
for hadronic interactions (jp interactions) becomes larger 
than the cross section for pair production. The crossover 
energy between these two processes is between 10 and 
10 20 eV and is found to be more or less independent of 
the material. This may be explained by the fact that 
the increase in a lp /<j ee for heavier nuclei is canceled out 
by the decrease of Elpm as X drops, as underlined in 
[16, 17]. This crossover energy is highly dependent on the 
model used to calculate the photonuclear cross section. 

This cross section rise reduces the mean free path of 
photons in the electromagnetic shower and therefore re- 
duces the total shower length. In addition, a hadronic 
and muonic component appears in the shower which will 
no longer be purely electromagnetic. Part of the energy 
of the electromagnetic shower leaks to this hadronic com- 
ponent and this will further decrease the shower length. 

Our goal here is to study the shower length and en- 
ergy under these two processes via a complete and re- 
alistic Monte Carlo simulation of underground electron- 
induced showers. One can mention that the need for 
detailed simulation of photonuclear interactions was al- 
ready mentioned in [22]. 

In our study we used TIERRAS [23] to simulate under- 
ground showers. TIERRAS is an adaptation of AIRES 
(Air-shower Extended Simulations) [24] for material like 
rock or ice. This simulation code provides a complete and 
realistic Monte Carlo treatment of the interactions which 
takes place within the shower development in dense mat- 
ter. It implements Migdal's theory of the LPM effect 
exactly and solves recursively the Migdal's cross section 
parameters [14]. 

The definition of standard rock used in this study refers 
to rock with constant density p = 2.65 g.cm -3 (density of 
continental earth crust), effective Z = 11 and Zj A = 0.5 
(i.e., isoscalar nuclei) [23]. The density of the earth crust 
used by the HiRes collaboration in [7] is 2.80 g.cm -3 . 
We do not expect this small discrepancy to impact our 
conclusions because the cross over energy between PN 
and LPM cross sections is essentially independent of the 
density. 



3 




FIG. 1. Pair production under LPM regime (LPM, black) and 
photonuclear interaction (PN, red) cross sections in standard 
rock used in TIERRAS [23]. 



The photonuclear interactions are processed with 
QGSJETII [26] and the cross section is given by a 
paramctrization and extrapolation of experimental data 
given in [27]. One should note that QGSJETII does not 
take into account the production of charmed mesons and 
associated prompt muons production. 

Cross sections used in this study are shown in Fig.l 

As shown in [16], another effect plays a role in the 
underground shower development: the electronuclear in- 
teraction. This effect is currently not taken into account 
by TIERRAS but it would be interesting to implement 
it in the future. 



III. MONTE CARLO STUDY 
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FIG. 2. Longitudinal Profile for shower induced by a 100 EeV 
electron with "LPM only" in the top and with "LPM and pho- 
tonuclear" processes in the bottom. The all particle profile 
is shown in red. The different components, namely electrons, 
photons, pions and muons are respectively displayed in blue, 
green, black and yellow. 



Using TIERRAS, we simulated underground showers 
induced by electrons with an energy between 10 17 eV 
and 10 20 eV (which is the range where the transition be- 
tween the LPM and PN cross sections is expected). In 
these simulations, we propagate the showers first using 
standard brcmsstrahlung and pair production processes 
(later called Bcthe- Heftier, "BH case"), then accounting 
for the LPM effect ("LPM only" case) and finally en- 
abling the photonuclear effect ("LPM+PN case"). 

The hadronic component produced by photonuclear in- 
teractions is also simulated even though we do not expect 
it to propagate over large distances: we thus propagate 
and follow electrons, photons, muons and hadrons. 

As in the case of air-shower simulations, for UHE pri- 
maries, the number of particles that are produced is large 
enough to make it impossible to propagate all the secon- 
daries. The simulations are made possible thanks to a 
sampling algorithm: the Hillas thinning algorithm [25]. 
One of the drawback of this algorithm is to introduce 
artificial fluctuations in the distribution of different ob- 



servables such as the longitudinal profile, the lateral dis- 
tribution, or the energy of the shower. In order to get 
smoother distributions we used a low thinning factor of 
10~ 6 (given by E t h /E pr i m , where E t h is the energy below 
the thinning algorithm is applied and E pi i m the energy 
of the primary) . In order to control the fluctuations in- 
duced by sampling, we also set the statistical weighting 
factor Wf to 0.2 [24]. For more information about the 
impact of the thinning algorithm on shower simulations, 
see [24]. 

On Fig. 2 one can see the longitudinal profile for 
10 20 eV electron showers, averaged over 100 showers, for 
each type of particle, and for the following cases: with 
LPM suppression only ("LPM only") and with LPM and 
photonuclear effects (" LPM+PN" ) . One can see that the 
"LPM only" showers will develop along much larger dis- 
tances than the "LPM+PN" showers. One can also see 
the apparition of hadrons and muons in the shower in the 
"LPM+PN" case highlights the fact that a part of the 
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FIG. 3. Energy spectrum of muons at the mean maximum 
development of the underground shower 



total shower energy is transmitted to hadronic subshow- 
ers and that a muonic component builds up and propa- 
gates much further away than the rest of the shower and 
carries a non negligible part of the total shower energy 
This is of course of interest for an experiment sensitive 
to muons such as the Pierre Auger Observatory. Figure3 
shows the energy spectrum of muons above 1 GeV at 
the maximum of development of the shower. Taking into 
account charmed meson production producing prompt 
muons would harden this spectrum. 

If one would look at the longitudinal profiles of in- 
dividual showers one would see important fluctuations 
from shower to shower. This is a characteristic behavior 
of electromagnetic showers in the LPM regime, as under- 
lined in [17] and indeed confirmed by our simulation. As 
we can see on Fig. 4, there is an important difference be- 
tween the mean longitudinal profile and its median. This 
is due to the fact that the mean is dominated by a few ex- 
treme values coming from a few showers. If one looks at 
the first and third quartile, we can see that the interquar- 
tile range is large: ~ 1 decade at the shower maximum 
development for the LPM case and up to ~ 2 decades 
for the LPM+PN case where the fluctuations are more 
important close to the maximum. This shows that one 
should be careful with the interpretations of the mean 
longitudinal profiles. Afterward, in order to avoid bias 
by averaging on all the showers, our calculations of the 
probability of an underground shower emerging in the 
atmosphere are made on a shower to shower basis. 

In order to estimate the shower length we need to de- 
fine an estimator of the latter. We adopted the definition 
of the "shower length" as given by Alvarez-Muniz and Zas 
in [28] where it is defined as the length over which the 
number of particles in the shower is larger than an arbi- 
trary fraction of the maximum number of particles in the 
shower. We choose here to set this fraction to 1/10. Re- 
sults arc shown in Fig. 5, where we plot the shower length 
as a function of the energy of the primary electron. 
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FIG. 4. Mean longitudinal profile (black), median (red), 
first and third quartile (blue) for shower induced by a 100 EeV 
electron with "LPM only" in the top and with "LPM and pho- 
tomiclear" processes in the bottom. 



One can notice that photonuclear processes make no 
difference in the shower length below 10 18 eV because the 
photonuclear cross section is dominated by radiative pro- 
cesses at these energies as can be seen from the cross sec- 
tions (Fig.l). The shower length increases strongly with 
energy when the LPM effect is applied without photonu- 
clear effect. This increase is severely reduced when both 
LPM and photonuclear effects arc taken into account. 

These results arc in good agreement with [16] in which 
the same behavior is predicted from a phenomenological 
approach. 

In a second study, we then consider the probability 
that the EM energy of the emerging shower tail exceeds 
a given detection threshold. A reasonable estimate of 
what would be the energy threshold for an array of sur- 
face detectors such as the Pierre Auger SD for detecting 
quasi-horizontal grazing showers emerging inside the ar- 
ray boundaries should be close to the 10 17 eV. We thus 
fix the energy detection threshold to this value, although 
more realistic results would require a complete detector 
simulation. 
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FIG. 5. Shower length versus primary electron energy. In- 
cluding "Bethe-Heitler" (squares), "LPM only" (triangles), 
"LPM+PN" (dots). 



We used the same set of simulated showers as above 
(i.e. showers from primary electrons with energy between 

10 17 and 10 20 eV). 

On Figure 6 we represent the probability that an un- 
derground shower with primary energy of 10 18 eV and 
10 20 eV emerges with a remaining energy above 10 17 eV 
after crossing a given amount of rock. 

One can see that the probability for a shower emerg- 
ing with an energy above 10 17 eV is higher in the "LPM 
only" case than in the two other cases. Photonuclear 
processes increase the energy loss rate in a non negligible 
way. If we compare the "LPM only" and "LPM+PN" 
case one can notice that at 10 20 eV the total probabil- 
ity is up to 2 times higher in the LPM case while at 

10 18 eV this total probability is only ~ 1.2 times higher. 
At 10 18 eV the overestimation of the emerging probabil- 
ity is not really important, but above 10 19 eV, PN in- 
teractions start to make a difference between the "LPM 
only" and "LPM+PN" case. 

In [7], the High Resolution Fly's Eye collaboration 
published a limit for up-going v e inducing LPM show- 
ers in crust. In their article, the authors define et as the 
transmission probability, which is the probability for a 
shower to emerge into the atmosphere with more than 10 7 
particles (above a low energy threshold). This probabil- 
ity can easily be extracted along with the shower length 
in our simulations. The energy dependence we obtain for 
et shows the same trends as the curves corresponding to 
the mean shower length. This probability is equivalent 
to the ratio between the mean shower length and the to- 
tal thickness of rock considered (by defining the shower 
length as the length where the shower has more than 10 
particles). Here we are considering a total rock thickness 
of 600 m and since we do not know the energy threshold 
used by the HiRes collaboration, we choose to fix this 
threshold at 500 MeV regardless of the type of the par- 
ticle, assuming that particles with lower energy will not 
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FIG. 6. Probability of an underground shower with primary 
energy of 10 18 eV (in the top) and 10 20 eV (in the bottom) to 
emerge with an energy above 10 17 eV after crossing a given 
amount of rock (in meters on this graph) for the three cases: 
"BH", "LPM only" and "LPM+PN". 



contribute significantly to the detection. 

By comparing the values of e t in the "LPM only" case 
and in the "LPM+PN" case, one can say that neglecting 
PN interaction at 10 18 eV only leads to an overestimation 
of e t of ~ 10%, at 10 19 eV this overestimation grows to 
- 40% and at 10 20 eV to - 110%. 

One concludes that neglecting the photonuclear effect 
causes an overestimation of the emerging shower energy 
and of the emerging number of particles above 10 eV. 
The detection probability is related to e*, which is itself 
related to the total energy and number of particles emerg- 
ing into the atmosphere. Overestimating the emerging 
energy as well as the shower length will lead to an arti- 
ficial enhancement of e* above 10 19 eV and therefore of 
the experimental sensitivity to this type of event. 
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FIG. 7. Probability of an underground shower to emerge with 
a number of particle above 10 7 versus primary energy, normal- 
ized to 600 m of rock. "LPM only" in blue and "LPM+PN" 
in red. 

IV. CONCLUSION 

Accounting for photonuclear interactions into our sim- 
ulations, we come to the conclusion that the LPM effect 
still enhances the shower length in rock but this enhance- 
ment is limited in a non negligible way. The underground 
shower energy loss is also impacted by the photonuclear 
interactions, as we can see in the decrease of the prob- 
ability to emerge above a given energy. In counterpart 
a non negligible muonic component will propagate over 



large distances under rock carrying a part of the shower 
energy. It appears that one cannot neglect the photonu- 
clear effect in simulations of underground showers where 
its cross section dominates the radiative processes cross 
sections above 10 19 eV. 

The development length of v e induced electromagnetic 
showers, the emerging shower energy and therefore the 
exposure to up-going v e event have been overestimated 
in [7] (and more particularly above 10 19 eV. The corre- 
sponding v e flux limit should thus be reevaluated by the 
authors by taking photonuclear processes into account. 

Another important fact may impact further the devel- 
opment of LPM showers. As outlined by some authors 
[16, 29], the radiated bremsstrahlung photon may inter- 
act via the photonuclear effect before being actually fully 
formed and this could further limit the the LPM inter- 
action length. We have not yet implemented this effect 
into our simulation and it is so far not taken into account 
in TIERRAS, AIRES and CORSIKA [30] codes. Taking 
into account this effect would further reduce the shower 
length and hence increase the discrepancy between LPM 
case and a treatment including all processes. 

Further work is hence needed to account in full detail 
for emerging showers and evaluate the sensitivity of UHE 
cosmic ray observatories to the flux of v e searching for 
these intriguing events. 
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